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Abstract 

Through  collaboration  with  Dr.  K.  H.  Schoenbach  of  Texas  Tech 
University  the  plasma  focus  opening  switch  (PFOS)  was  revised  to  answer 
basic  questions  as  to  the  utility  of  the  concept.  To  estimate  the  plasma 
temperature  and  classical  resistivity  a  soft  X-ray  spectrometer  and  X-ray 
pinhole  camera  were  developed.  The  temperature  was  estimated  from  a  coronal 
model  to  range  between  0.4  to  0.5  keV  for  either  a  nitrogen  or  neon  impurity 
(1  to  2%)  in  deuterium  at  3  torr.  Strong  pinches  were  observed  in  pure  neon 
(0.6  torr)  with  an  electron  temperature  in  the  same  range.  The 
corresponding  classical  resistance  of  the  pinch  is  9mfi  whereas  500mfi  is  more 
consistent  with  output  voltage  pulse  and  current  flow  at  interruption 
indicating  anomalous  resistivity  is  present.  A  one-dimensional  two-fluid 
computer  code  has  been  developed  to  model  anomalous  resistivity  in  the  pinch 
phase  and  preliminary  results  are  consistent  with  the  snowplow  model.  The 
final  analysis  of  the  plasma  focus  particle  beam  generation  experiments  was 
completed  and  a  strong  correlation  was  found  for  the  beam- target  model  as 
the  mechanism  for  neutron  production  in  the  Illinois  plasma  focus  device. 


As  indicated  on  the  milestone  chart  for  the  1981-82  proposal  (Figure  1) 
the  major  research  tasks  were  the  development  of 

1)  Soft  X-ray  diagnostics 

2)  Laser-interferometer  for  density  measurements 

3)  Computer  modeling  of  the  plasma  focus 

4)  High-Z  gas  operation  of  the  plasma  focus. 

Due  to  budget  limitations  only  the  first  item  was  completed  and  the  rest  are 
in  various  stages  of  development.  Even  to  achieve  this  stage,  equipment  and 
computer  time  contributed  by  other  funding  sources  had  to  be  utilized.  Also 
collaboration  with  Prof.  K.  Schoenbach  from  Texas  Tech  University  (TTU) 
caused  the  direction  of  the  research  to  be  modified.  This  change  emphasized 
the  measurement  of  the  plasma  resistance  for  which  the  soft  X-ray 
diagnostics  are  the  most  crucial.  Thus  all  things  considered  good  progress 
was  made  which  will  pave  the  way  for  completing  the  revised  goals  (as  shown 
in  Figure  2  for  FY  1982-83)  in  a  timely  manner, 
c)  Status  of  Plasma  Focus  Opening  Switch  Research 
1)  Overview:  collaboration  with  Prof.  K.  H.  Schoenbach  of  Texas  Tech 

University 

After  meeting  Prof.  Schoenbach  at  the  3rd  Pulsed  Power  Conference  in 
Albuquerque,  NM  in  June  1981  and  after  having  had  several  telephone  conver¬ 
sations  with  him  since  that  meeting  it  was  realized  that  the  possibility  of 
fruitful  collaboration  existed.  This  was  strongly  encouraged  by  the  grant 
monitor  at  that  time,  Lt.  Col.  A.  K.  Hyder,  so  in  January  1982,  G.  Gerdin 
(principal  investigator  on  this  grant)  traveled  to  Texas  Tech  University  for 
a  three  day  visit  to  discuss  the  possible  nature  of  this  collaboration. 
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Prof.  Schoenbach  is  an  ideal  man  with  whom  to  collaborate  since  he  had 
been  involved  with  dense  plasma  focus  (DPF)  research  for  about  ten  years  and 
is  the  author  of  over  ten  papers  on  the  subject.  Presently  he  is  directing 
a  project  on  diffuse  discharge  opening  switches  and  so  has  experience  in 
both  DPF  and  opening  switch  physics. 

A  previous  unsuccessful  attempt  to  use  the  plasma  focus  as  an  opening 
switch*  was  regarded  as  not  a  true  test  of  a  plasma  focus  as  an  opening 
switch  (see  Appendix  A).  That  is,  it  did  not  appear  the  device  had  been 
operated  in  an  optimum  manner  since  no  capacitive  discharge  data  showing 
strong  interruptions  were  presented^  and  the  initial  I0  was  very  low.^ 
Furthermore  the  circuit  used^  represented  a  brute  force  approach  and  the 
results  indicated  a  more  sophisticated  circuit  or  gas  handling  technique 
would  be  required  (see  Appendix  A).  So  a  research  program  was  devised  to 
provide  experimental  tests  to  gain  an  understanding  of  the  DPF  interruption 
physics  under  ideal  conditions  of  a  capacitive  energy  drive  and  static  gas 
fill. 

The  research  program  will  answer  these  four  basic  questions: 

O 

a)  Is  the  interruption  (Figure  3)  inductive  (L)  or  resistive  in 
nature? 

b)  Can  the  power  be  tapped  into  a  load? 

c)  Can  the  timing  of  the  interruption  be  controlled? 

d)  Can  the  operation  be  made  repetitive? 

The  resulting  program  is  discussed  in  detail  in  Appendix  B  and  presented  in 
Figure  2. 


The  work  performed  during  1981-82  will  have  the  most  impact  on  the 
first  year  of  this  plan.  The  main  goal  of  the  first  year  of  this  program  is 
to  measure  the  resistance  of  the  plasma  focus  and  compare  it  with  classical 
resistivity.  To  determine  the  latter,  experimental  measurements  of  the 
electron  temperature,  the  plasma  charge  state  and  shape  are  the  most 
crucial. ^  These  parameters  can  be  resolved  by  soft  x-ray  diagnostics  and 
the  nature  and  results  of  these  measurements  will  be  discussed  in  section 
C2. 

The  computational  effort  has  centered  about  the  development  of  1-D 
two  fluid  MHD  computer  code  to  model  the  collapse  (interruption)  phase  of 
the  DPF  and  the  status  of  these  calculations  will  be  discussed  in  section 
C3. 

The  final  analysis  of  the  beams  generated  by  the  plasma  focus  based  on 
results  reported  last  year  will  be  presented  in  section  C4. 


C2:  Soft  X-ray  Analysis  of  a  Dense  Plasma  Focus 
General  Introduction  and  Purpose 

An  experiment  to  detect  and  analyze  soft  x-ray  radiation  coming  from 
the  pinch  of  the  Illinois  Dense  Plasma  Focus  (DPF)  has  been  designed, 
constructed  and  performed  during  the  summer  and  fall  semester  1982. 

The  experiment  was  designed  to  determine  the  electron  temperature  of  the 
D.P.F.,  an  important  parameter  in  most  plasma  calculations.  The  detection 
of  the  x  rays  is  performed  by  silicon  solid  state  devices  (pin  diodes) 
conveniently  filtered  by  matched  metal  foils  in  an  arrangement  known  as 
Ross-Filter  Technique^-?.  The  analysis  of  the  experimental  data  is 
handled  through  a  computer  model  of  the  observed  plasma  to  give  the  electron 
temperature  as  a  function  of  an  experimentally  measured  parameter:  the 
ratio  between  the  detected  intensities  of  selected  spectral  bands6. 

Hardware 

The  x-ray  spectrometer  (XRS)  consists  of  a  long  (3.0  m)  evacuated 
aiming  tube  with  a  tightly  packed  array  of  three  pin-diodes  at  one  end.  An 
aluminized  Mylar  vacuum  window  at  the  other  end,  housed  in  a  revolving  ball 
valve,  couples  the  XRS  with  the  focus,  allowing  easy  aiming  of  the 
spectrometer.  The  interface  of  the  XRS  with  the  DPF  chamber  is  completed  by 
an  adapter  flange  which  permits  vertical  movement  of  the  spectrometer  and  by 
a  bellows  connection  which  allows  horizontal  movement.  An  extension  pipe 
protrudes  into  the  DPF  chamber,  as  close  as  2  inches  away  from  the  axis, 
allowing  excellent  space  resolution.  (Figure  4).  The  device  is  aimed 
through  the  combination  of  bellows,  adapter  flange  and  extension  pipe  and 
this  operation  is  performed  routinely  by  members  of  the  DPF  experimental 
group  before  each  data  run. 
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The  spectrometer  is  kept  under  vacuum  (<5y)  by  a  sieve-trapped 
mechanical  pump.  This  prevents  the  soft  x  rays  from  being  significantly 
attenuated  on  their  way  to  the  detectors.  The  spectrometer  is  also 
electrically  insulated  from  the  OPF  chamber  (which  rises  to  high  potential 
during  the  capacitor  bank  discharge  at  the  firing  of  the  focus). 

The  solid  state  detectors  used  in  the  experiment  are  Quantrad 
100-PIN-250  doubly-diffused  silicon  pin  diodes  (Figure  5).  They  were  chosen 
because  of  their  high  sensitivity  in  the  region  of  interest  (0.1  to  10  keV) 
and  fast  rise  time  (~2  nsec.)  which  allow  good  resolution  in  both  space  and 
time^. 

Each  pin  diode  is  filtered  by  a  different  metal  foil.  The  effective 
thicknesses  of  the  foils  (thickness/sine  of  angle  of  inclination)  are 
matched  so  that  their  attenuation  power  is  the  same  for  all  three  filters 
except  in  the  regions  between  absorption  k  edges.  (Figures  6,7,8).  The 
metal  foils  used  are 

1.  Beryllium  (1. l*10"lg/cm2;  k  edge  at  0.11  keV) 

2.  Aluminum  (4.1*10~3g/cm2;  k  edge  at  ~  1.5  keV) 

3.  Nickel  (4.9’10“4g/cm2;  k  edge  at  ~  8  keV) 

By  subtracting  the  signal  of  pin  diode  one  from  that  of  two  (2-1)  and  the 
signal  of  pin  diode  two  from  that  of  three  (3-2)  it  is  possible  to  obtain 
two  values  proportional  to  the  intensity  of  the  radiation  emitted  in  the 
regions  between  k  edges  (0.1  to  1.5  keV,  1.5  to  8  keV).  Pictures  of  the 
pin-diode  signals  are  included.  (Figures  9,  10). 


This  information  is  the  experimental  data  that  is  compared  to  the  curve 
generated  by  a  computer  simulation,  according  to  a  given  theoretical  model, 
to  yield  a  value  for  the  electron  temperature.  Future  improvements  of  tne 
detection  system  include  the  addition  of  two  more  semiconductor  detectors, 
for  a  total  of  five  pin  diodes,  thus  defining  four  adjacent  radiation 
bands. 

In  addition  to  the  pin  diode  XRS,  a  pinhole  camera  was  constructed  and 
used  to  detect  time  integrated  x-ray  radiation  from  the  focus.  (Figure  4). 

The  pinhole  camera  is  connected  to  the  chamber  port  opposite  to  the 
x-ray  spectrometer.  A  1  mil  thick  beryllium  window  is  used,  with  a  20  mil 
diameter  tantalum  pinhole.  The  camera  is  flushed  with  helium  during 
operation,  to  minimize  attenuation  of  x  rays.  No  x-ray  radiation  could  be 
detected  without  helium  flushing,  thus  indicating  the  predominance  of  x  rays 
from  the  soft  region  in  the  focus  emission  and/or  the  lack  of  film 
sensitivity  to  harder  x  rays. 

Software 

The  computer  simulation  of  the  DPF  plasma  x-ray  emission  uses  a  steady 
state  optically  thin  corona  model  approximation.  This  is  consistent  with 
previous  studies  on  the  subject. 8,9  code  has  been  developed  to 
calculate  the  x-ray  emission  from  the  focus  and  the  frequency-integrated 
response  of  the  filtered  pin  diodes  as  a  function  of  electron  temper¬ 
ature.^  The  following  input  data  are  required: 

1.  Thickness  of  filter  foils  and  absorbing  power  as  a  function  of 

frequency  (energy). 

2.  Composition  of  the  pin  diodes. 
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3.  Cross  sections  for  collisional  ionization  and  radiative  recombina¬ 
tion 

4.  Probabilities  of  excitation  and  de-excitation  for  various  states.  ^ 

5.  Information  on  line  radiation  emitted  in  the  region  of  interest. 

The  code  evaluates: 

1.  Charge  state  equilibrium  conditions. 

2.  Continuum  radiation  (Bremstrahlung  +  recombination)^, 13 

3.  Line  radiation  (decay  of  excited  states). 

4.  Response  curve  of  filtered  pin  diodes. 

5.  Response  of  pin  diodes  to  the  x-ray  radiation  emitted  by  the 
plasma  defined  above. 

6.  Differences  in  magnitude  between  responses  (band  strengths)  and  ratios 
of  band  strengths. 

The  final  output  is  a  table  of  ratios  of  band  strengths  vs.  electron 
temperature,  to  be  compared  with  the  experimentally  observed  ratios. 

A  steady  state  corona  model  is  presently  used  in  the  code. 5, 8  ^ 
more  sophisticated  version  of  the  corona  model  is  now  being  considered  for 
implementation,  which  uses  more  recent  and  accurate  values  of  cross 
sections. 6,14-16  Also  a  local  thermodynamical  equilibrium  (LTE)  model 
is  going  to  be  used  to  evaluate  what  differences,  if  any,  the  LTE  assumption 
introduces  in  the  determination  of  the  plasma  state  and  the  electron 
temperature.  Eventually  a  full  radiation-coil isional  time  dependent  code 
will  be  used  in  conjunction  with  an  MHO  code  to  model  the  x-ray  emission 
during  the  pinch  phase  of  the  plasma  focus. 16-18 
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Results 

Preliminary  results  based  on  simple  steady  state  corona  model  show  good 
consistency,  yielding  temperature  values  of  about  0.4  to  0.5  keV,  regardless 
of  the  type  of  impurity  gas  present  (neon  vs.  nitrogen).  (Figure  11).  This 
result  is  encouraging  since  most  of  the  detected  x-ray  emission  from 
nitrogen  seeded  plasma  is  of  the  continuum  type,  as  opposed  to  a  neon 
contaminated  plasma,  where  the  neon  line  radiation  dominates  the  detected 
emission.  Roughly  the  same  temperatures  are  expected  in  both  cases,  given 
the  small  amount  of  impurity  gas  present  (about  2%).  (Figure  11).  Strong 
pinches  were  observed  with  pure  neon  as  the  working  gas  (0.6  torr)  and  with 
temperatures  about  the  same  as  reported  above. 

Copper  line  radiation  from  the  center  electrode  (hollow)  has  been 
determined  to  be  negligible,  on  account  of  the  following  results: 

1.  X-ray  emission  increases  by  a  factor  greater  than  10  when  impurities 
(neon  or  nitrogen)  are  present.  Almost  no  signal  is  detected  by  the 
filtered  pin  diodes  for  pure  deuterium  fills. 

2.  X-ray  pinhole  camera  pictures  show  x-ray  emission  mostly  from  the  pinch 
region.  (Figure  12).  Compare  with  Ref.  4  (solid  center  electrode). 

The  time  integrated  pictures  of  the  x-ray  emission  from  the  focus  also 

indicates  an  upper  limit  of  1-2  mm  for  the  diameter  of  the  pinch.  When  the 
absolute  response  of  the  pin  diodes  (supplied  by  the  manufacturer)  is  used, 
together  with  the  nominal  solid  angle  of  view,  the  measured  size  of  the 
pinch,  and  the  actual  signal  observed,  the  electron  density  can  be 
estimated.  The  values  of  the  electron  density  calculated  in  this  way  are 
between  6*10^  and  3*10^  electrons/cm^  for  an  optically  thin  plasma. 


A  simple  check  of  this  estimate  can  be  made  assuming  magnetic  to  thermal 

pressure  balance  (0=1)  and  using  experimentally  evaluated  parameters  (pinch 

current  =  300  kA,  pinch  radius  =  0.5-1  mm,  electron  temperature  =  500  eV). 

This  yields  values  for  the  electron  density  in  the  same  range  (8-10^ 

to  5- 10^9  electrons/cm^)  depending  on  assumptions  for  the  ion  temperature 

(T^=0  or  Tj=Te)  and  plasma  radius.  These  values  are  consistent  with 

19 

measurements  of  plasma  focus  densities  using  laser  interferometric  , 
and  laser  scattering*^  techniques  on  devices  of  similar  energy.  This 
fact  supports  the  validity  of  the  model  used  in  evaluating  the  experimental 
data  and  specifically  of  two  assumptions  that  have  been  made: 

1.  Line  radiation  from  other  impurities  (such  as  electrode  material)  is 
not  important  (the  model  ignores  it). 

2.  Self  absorption  of  radiation  by  the  plasma  is  negligible  (the  model 
uses  the  optically  thin  approximation). 

Also  it  is  possible  to  calculate  the  classical  (Spitzer)  resistance  of 
the  pinch  from  the  evaluated  dimension  of  the  x-ray  image  of  the  pinch. 

This  is  roughly  9  mft  at  a  temperature  of  500  eV.  This  value  is  about  fifty 
times  lower  than  the  observed  0.5ft  value  (from  current  interruption 
measurements)  suggesting  the  possibility  of  non-classical  (anomalous) 
resistivity  playing  an  important  role  in  the  focus  current  interruption  at 
pinch  time. 

Conclusions 

A  combination  x-ray  spectrometer  and  pinhole  camera  has  been  built  and 
tested.  X-ray  spectral  analysis  of  radiation  from  the  Illinois  dense  plasma 
focus  has  been  performed  to  determine  the  electron  temperature  of  the 
pinching  plasma.  Preliminary  results  indicate  a  temperature  in  the  range 


of  400  to  500  eV.  These  results  are  supported  by  cross  checks  on  electron 
density  calculations.  Improvements  have  been  planned  and  are  being 
implemented  in  both  the  experimental  hardware  and  software. 

C3;  Theoretical  and  Computational  Effort 

Although  the  research  is  primarily  an  experimental  study,  about  one 
third  of  the  effort  is  directed  toward  developing  models  for  the  plasma  focus 
opening  switch  behavior  in  circuits  with  realistic  loads  to  obtain  an 
understanding  of  the  coupling  of  the  circuit  and  plasma.  Specifically  the 
role  of  fluctuations  in  creating  anomalous  resistivity  at  pinch  time  appear 
particularly  relevant.  The  results  of  the  development  of  theoretical  models 
will  be  compared  with  experimental  results  on  a  frequent  basis  to  suggest 
further  experimental  tests  and  new  ways  toward  achieving  a  repetitive  opening 
switch. 

Two  computer  codes  have  been  developed.  The  first  is  a  zero-dimensional 
circuit-dynamical  model  where  the  macroscopic  radial  motion  of  the  plasma  is 
predicted  self  consistently  using  the  snow  plow  model,  the  Leontovich- 
Osovets'  equation,  and  the  circuit  equation. 21  classical  resistivity^ 
has  been  included  in  the  model  and  the  results  will  be  compared  with  measured 
circuit  parameters  to  try  to  determine  the  role  of  sheath  motion  ( L * ) .  This 
code  has  been  used  to  predict  the  magnitude  of  the  voltage  pulses.  These  are 
arising  from  sheath  motion  and  they  are  found  to  increase  with  magnitude  of 
current  flow  at  interruption^  in  the  same  manner  and  magnitude  as 
observed  experimentally. 24 

To  do  a  more  detailed  check  on  the  theoretical  possibility  of  anomalous 
resistivity  in  a  plasma  like  the  plasma  focus,  a  one  dimensional  (cylindrical 
geometry)  computer  model  (DPFR)  of  the  DPF  has  been  written^  in  order  to 
compare  experimental  results  with  theoretical  anomalous  transport 


coefficients  for  various  microinstabilities. 25  The  equations  are 
integrated  using  the  Lax-Wendroff-Richtmeyer  finite  difference  scheme2^ 
with  the  help  of  flux -corrected  transport. 27  This  approach  has  been 
successfully  used  in  modeling  theta  pinch  experiments23  has  been  adapted 
to  DPF  geometry.  The  code  calculates  the  density,  temperature  and  axial 
electric  field  in  the  DPF  as  a  function  of  time  and  includes  the  circuit 
equation  self  consistently.  (Figure  13). 

Previous  DPF  models29»30  have  been  based  on  pure  MHD  equations  which 
cannot  include  the  effects  of  turbulence,  such  as  anomalous  resistivity,  and 
therefore,  are  unable  to  accurately  predict  the  characteristics  of  DPF 
impedance. 

For  an  initial  pinch  current  of  450  kA  with  the  ion  acoustic  instability 
included,  the  code  predicts  a  peak  plasma  density  of  3.0  x  10liJ  cm-3, 
peak  electron  temperature  of  1.9  keV,  dl^g/dt  =  -2.2  x  1012  amps/sec 
(where  I^g  is  the  main  bank  current)  and  a  peak  resistivity  of  about  1000 
X  Spitzer.  (The  experimentally  observed  value  of  dl^g/dt  during  the 
radial  collapse  is  about  -  1.9  x  10^2  amp/sec;  see  Figure  3.)  Calculated 
profiles  of  the  density,  electron  temperature,  ion  temperature,  and  the 
magnetic  field  are  shown  in  Figure  14  at  about  one  nanosecond  before  the 
electrons  run  away. 

The  code  predicts  the  collapse  time,  tc,  (at  constant  fill  pressure) 
to  be  inversely  proportional  to  the  initial  current,  I0,  (see  Figure  15)  as 
consistent  with  the  snowplow  model^l.  The  short  collapse  times  are  due 
to  the  initial  density  distributions  (Figure  16)  which  were  selected  to  match 
boundary  conditions  and  shorten  the  computation  time.  More  realistic  initial 
density  profiles  will  be  run  in  the  near  future.  The  high  ratio  of  Te/T^ 
shown  in  Figure  14  are  due  to  the  lack  of  the  effects  of  viscosity  in  the  code 
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at  present.  The  addition  of  classical^!  viscosity  into  the  code  is  part 
of  the  present  effort.  Thus  the  computational  effort  is  off  to  a  good 
start. 

C4:  Summary  of  Particle  Beam  Scaling  Measurements 

The  results  reported  in  the  last  annual  report^, 33  have 
implications  on  two  aspects  of  DPF  phenomena: 

1.  the  neutron  production  mechanism,  and 

2.  the  acceleration  mechanism. 

These  aspects  will  be  treated  in  this  section. 

Since  the  total  number  of  particles  accelerated  and  their  energy  spectra 
have  been  determined  one  can  use  this  information  to  predict  the  magnitude 
and  the  scaling  of  the  axial  beam  target  neutron  yield  with  IMB  if  a  few 
simple  assumptions  are  made  about  the  target.  If  the  target  is  assumed  to  be 
the  cold  gas  ahead  of  the  beam  then  it  is  independent  of  bank  parameters. 
Similarly,  if  one  assumes  that  a  hot  target  (plasma)  of  density  and  length 
independent  of  I^g,  which  is  roughly  consistent  with  several 
experiments. 19.20,34-36  ^ven  slowing  down  of  the  deuterons  could  be 
neglected  for  Te  >  500  eV  for  targets  less  than  10  cm  in  length.  Thus  in 
calculating  the  axial  beam  target  neutron  yield  the  following  assumptions 
have  been  made: 

1.  The  ion  beam  energy  spectrum  has  the  same  power  law  dependence  as 
the  electron  beam  energy  spectrum  for  a  given  I^g. 

2.  N+/Ng  =  (M+/Me)1/2. 

3.  The  ion  and  electron  energy  beam  spectra  are  the  same  at  all  angles 
to  the  axis  where  the  beams  exist. 

4.  The  target  parameters  are  independent  of  I^g. 
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5.  The  slowing  down  in  the  cold  gas  target  is  that  for  deuterons  in  3  torr 
D2  gas.37 

6.  The  slowing  down  in  the  hot  plasma  target  is  negligible. 

The  predictions  of  this  model  are  compared  with  the  experimental  data  in 


Table  I. 


Table  I 


Phenomena 


Y  *  x 
n  XMB  *  A 

Target  thickness 


Yn  at  12.5  kJ 

Anisotropy 

(Yn0o/Yn90o) 

Neutron  pulse 
width  at  3.2  meters 


Experiment 


4.6  +  0.3 


Cold  Target  Hot  Target 


5.2  +  1.0 


4.8  +  1.0 


15  cm  @  3  torr  2.9  x  lOiadeuteron 


3.8  x  10* 
1.7  +  0.2 


90  ns 


2.4  x  10- 


+60  ns 


3.1  x  10" 


+30  ns 


+Assumes  axial  beam  motion  only  and  includes  20  ns  pulse  widtn. 

The  scaling  and  magnitude  of  the  beam  target  neutron  yield  is  consistent 
with  either  target  model  (Figure  17)  and  the  electron  beam  scaling  parameters 
reported  previously^SS.  However,  neither  target  model  is  adequate  to 
explain  the  neutron  pulse  width.  The  long  pulse  width  oDserved  experimentally 
could  be  due  to  trapping  of  part  of  the  beam  in  some  magnetic  structure  of  the 
focus  as  evidenced  by  the  high  ratio  of  DT  neutrons  to  DD  neutrons  observed  in 
other  DPF  experiments. 34»38  Description  of  such  trapping  is  beyond  the 
scope  of  these  target  models.  The  observed  anisotropy  is  consistent  with  only 
the  hot  target  model. 


Due  to  the  assumptions  made  in  the  calculations,  the  agreement  in  the 
magnitude  and  scaling  of  the  beam  target  model  prediction  of  neutron  yield 
and  the  experimental  result  only  means  that  the  axial  beam  target  model 
cannot  be  ruled  out.  Probably  the  most  questionable  assumption  is  that  the 
deuteron  beam  energy  spectrum  is  the  same  at  all  angles  (to  15*)  with  respect 
to  the  device  axis.  This  is  apparently  true  for  the  electron  beam  which  is 
not  always  on  axis  (presumably  due  to  the  frozen  hose  instability^)  and 
yet  always  gives  the  same  power  law  spectrum  (within  the  error).  The  large 
uniform  patterns  of  tracks  observed  in  the  SSNTD  converter  layer  technique4^ 
(which  has  a  lower  energy  threshhold  of  about  150  keV)  over  a  19*  angle  with 
respect  to  the  axis  also  supports  the  assumption.  And  finally,  the  deuteron 
energy  spectra  are  roughly  independent  of  angle  for  a  28  kJ  L)PF  at  5  torr 
02  in  the  energy  range  from  0.9  to  4.5  MeV.36. 

Acceleration  Mechanism 

It  is  clear  from  Figures  14  and  15  of  ref.  33  that  the  electron  (and 

O 

hence  ion  beam)  occur  at  the  peak  in  |  IMgj  indicating  that  current 
interruption  plays  a  role  in  the  acceleration  process.  However,  estimates  of 
the  induced  voltage  across  the  pinch  from  the  observed  voltage  surge  across 
the  parallel  plate  transmission  lines  between  the  capacitor  banks  and  the  DPF 

O 

device  and  LI  (where  L  is  estimated  from  pinhole  camera  pictures)  indicate 

O 

that  only  about  100  kV  appears  across  the  pinch  due  to  |  I^g|  phenomena. 

This  voltage  is  much  too  low  to  explain  the  high  particle  energies  observed 

O 

by  simple  diode  acceleration  if  current  interruption  (I)  was  the  only  source 
of  the  applied  potential. 
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It  has  been  suggested4*  that  motion  of  the  current  sheath  could 

O 

cause  induced  voltages  much  higher  than  LI  effects  and  the  computational 
results  of  Kondoh  and  Hirano2*  contain  estimates  of  the  magnitude  of  this 
effect.  The  system  of  equations  used  by  the  latter  authors  has  been 
programmed  and  used  to  estimate  the  induced  voltage  for  the  Illinois  DPF 
device.  The  final  radial  collapse  (r0  =  2.5  cm)  was  modeled  and  the 
initial  conditions  were  determined  experimental ly  by  placing  a  coaxial  short 
across  the  open  end  of  the  plasma  focus  electrode  and  measuring  the  circuit 
characteristics  when  the  device  was  evacuated.  At  3  torr  Dg  pressure,  the 
code  estimates  total  induced  voltages  range  up  to  ~ 100  kV  for  a  peak  device 
current  of  560  kA.  Hence  it  appears  that  even  sheath  motion  cannot  provide 
sufficient  potentials  to  explain  the  observed  particle  energies  (which 
apparently  range  above  1  MeV  in  this  device48)  by  the  simple  process  of 
diode  acceleration. 

Conclusions 

7 

The  energy  spectra  of  the  fast  ions  and  electrons  emitted  by  the  plasma 
focus  in  opposite  directions  (Figure  1  of  ref.  33)  are  observed  by  direct 
methods  to  have  the  same  power  law  energy  dependence  at  a  device  current  at 
pinch  time,  I^g,  of  560  kA.  The  scaling  of  the  absolute  value  of  the 
exponent  of  the  power  law  was  found  (for  electrons)  to  decrease  with 
increasing  Ip|g  (hardening  of  the  spectrum)  and  the  values  found  to  be 
consistent  with  the  electron  energy  spectra  inferred  from  hard  X-ray  spectral 
measurements  on  several  other  devices.42"46  it  is  interesting  that  power 
law  spectra  of  similar  exponents  have  also  been  observed  in  a  double  inverse 
pinch,46  cosmic  rays,47  and  solar  flares.48 
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The  magnitude  of  the  fast  electron  current  has  been  measured  by  a 
filtered  fast  Faraday  cup  and  found  to  scale  as  to  the  third  power. 

The  highest  primary  beam  current  observed  was  17  kA  for  a  bank  energy  of  12.5 
kJ.  Since  much  more  energetic  plasma  foci  exist  such  as  the  1MJ  device  at 
Frascati, it  would  be  very  interesting  if  similar  measurements  were 
performed  on  these  devices  to  see  if  the  scaling  observed  in  the  Illinois 
device  persists  at  high  bank  energy.  If  so,  the  hardening  of  the  energy 
spectra  indicated  by  these  measurements  and  others^-45  may  lead  to  new 
applications  of  the  DPF  as  a  pulsed  electron  beam  source. 

The  measured  beam  parameters  and  scaling  laws  can  be  used  to  predict  the 
scaling  and  magnitude  of  the  neutron  yield  due  to  a  beam  target  model.  If 
the  target  parameters  are  assumed  to  be  independent  of  I^b,  then  the 
agreement  between  beam  target  yield  and  the  observed  yield  is  quite  good 
with  respect  to  scaling  and  magnitude  but  the  predictions  of  the  model  do  not 
fit  the  observed  temporal  neutron  pulse  width.  It  is  possible  that  trapping 
of  all  or  part  of  the  deuteron  beam  by  magnetic  structures  in  the  plasma 
focus^4  could  explain  the  discrepancy. 

The  beams  are  observed  to  occur  at  the  peak  in  the  current  interruption 

O 

( Ifwie )  indicating  the  latter  has  a  strong  influence  on  the  acceleration 
process  but  the  magnitudes  of  the  voltages  generated  by  d/dt  (LI)  effects 
estimated  from  circuit  measurements  are  too  low  to  explain  the  acceleration 
process  by  simple  diode  acceleration. 

Since  the  beams  are  generated  during  the  power  pilse  at  the  location  of 
the  current  interruption  the  acceleration  mechanism  is  likely  to  be 
fundamentally  connected  to  the  mechanism  of  current  interruption.  From  the 
data  presented  here  plasma  diode  acceleration  formed  by  a  plasma  filled 
opening  switch  (current  interruption)  remains  a  strong  candidate  as  the 
accelerating  mechanism  and  hence  voltage  pulses  of  up  to  1MV  appear  to  be 
generated  in  the  switch. 
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Figure  1  Milestone  chart  for  the  1980-81  proposal. 


Figure  2  Milestone  Chart:  Plasma  Focus  Opening  Switch  Experiments 
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Figure  2  Milestone  Chart:  Plasma  Focus  Opening  Switch  Experiments 
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Figure  3 


Wave  form  of  the  device  current  of  the 
Illinois  plasma  focus  as  observed  with 
a  Rogowski  coil. 


■ 


Figure  4  Experimental  set-up,  x-ray  spectrometer  on  the  left,  dense  plasma  focus  chamber,  x-ray  pinhole 
camera  on  the  right.  Total  distance  between  x-ray  pin  diodes  and  focus  =  4.5  m,  diameter  of 
extension  pipe  =  1.0  cm;  sensitive  area  of  detectors  =  100  mm2. 
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Figure  5  Schematics  of  pin  diode  detector  circuit;  pin  diodes  are 
Quantrad  100  PIN  250;  bias  voltage  *  -300  volts;  a  =  angle 
of  Inclination  of  metal  foils. 
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Figure  8  Pin  diode  response  to  x-rays  with  nickel  filter  (4.9*10 

g/cnr),L-III  edqe  at  0.862  keV,  k  edge  at  8.33  keV  silicon 
edge  (pin  diode)  at  1.84  keV. 
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Figure  11  Ratio  of  band  strengths  vs.  temperature,  as  determined 
by  computer  simulation  (corona  model)  for  deuterium 
fills  with  2%  nitrogen  and  2%  neon. 
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Figure  15  The  dependence  of  radial  collapse  time  t  ,  on 

the  inverse  of  the  initial  pinch  currentc(l/I  ), 
as  computed  by  the  one-dimensional  two-fluid  0 


code.  The  snowplow  model  predicts  these  parameters 
(1/1.  O  to  be  linearly  related. 
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Figure  16  Initial  density  profile  as  used  at  present  in  the  two 
fluid  code. 
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Figure  17  Neutron  yield  in  the  Illinois  plasma  focus  and  the  predictions  of  various  models. 

4. 4+0. 3 

I)  Curve  a)  Is  the  least  square  fit  to  the  peak  neutron  yield  versus  I  (Y  “I  —  ). 

MB  n  MB 

II)  Curve  b)  is  the  least  square  fit  to  the  neutron  yield  for  the  shots  actually  used  to  deterlne  the 
electron-beam  generation  data.  The  triangles  are  the  data  with  representative  error  bars  shown  for 

4. 6+0. 3 

one  of  the  points.  For  these  data  Y  *  I  —  . 

n  Mb 


1 1 1 )Curve  c)  Is  the  result  of  the  beam  hot  target  model  for  a  target  thickness  of  2.9x10 


18  deuterons 


The  error  bars  are  those  associated  with  the  experimental  uncertainties  In  the  beam  parameters. 

4.8 

The  model  predicts  Y  «  I  . 

n  MB 


lv)  Curve  d)  Is  the  result  of  the  beam  cold  gas  target  for  a  target  of  3  torr  deuterium  gas  15  cm  long 
(target  thickness  2.9x1o'®  deuterons/cm^  as  In  III)  where  the  effects  of  slowing  down  on  the 
deuteron  beam  have  been  Included.  The  error  bars  are  as  In  ill)  and  this  latter  model  predicts 
„  „  ,  5.2 


Y  *  I  ' 
n  M 


A  previous  opening  switch  experiment1,  where  a  DPF  was  utilized  to 
interrupt  the  current  stored  in  a  large  inductive  coil  (see  Figure  Al),  was 
somewhat  unsuccessful  in  the  sense  that  the  interruptions  and  power  pulses 
generated  were  irregular  and  variable  in  magnitude. 

The  experiments  were  successful  in  that  several  focus  interruptions 
occurred  over  a  period  of  150  ysec  (a  factor  of  fifty  longer  than  normal) 
and  thus  these  results  show  that  a  long  series  of  power  pulses  may  be 
possible  if  only  the  moment  of  interruption  could  be  controlled. 
Understanding  the  nature  of  the  current  interruption  and  hence  possible 
means  of  achieving  firing  control  are  a  central  task  to  DPF -opening  switch 
research.  Some  progress  can  be  made  toward  this  goal  by  analyzing  the 
experiments  of  Ur.  Salge  in  more  detail. 

On  reviewing  the  results  of  Salge' s  experiments*  three  things  stand 

out: 

1)  The  voltage  across  the  device  is  always  a  few  kilovolts  during  the 
pulse. 

O 

2)  The  initial  time  derivative  of  the  device  current,  I0,  (as 
inferred  from  Figure  6  of  ref.  1)  is  low  (only  about  2*10*0 
A/s). 

3)  No  strong  interruption  occurs  at  the  initial  pinch  time  and  the 
maximum  interruption  is  only  20%  of  the  current. 

These  three  points  will  be  discussed  in  the  rest  of  the  appendix. 

The  first  point  illustrates  that  the  gas  in  the  Salge  DPF  experiments 
was  never  given  a  chance  to  recover  (recombine)  and  hence  each  successive 
pinch  and  interruption  depends  on  the  distribution  of  the  ionized  gas  left 
by  the  previous  one  (see  Figure  4  of  ref.  1).  Since  the  conventional  DPF 


even  under  ideal  conditions  is  not  very  reproducible  in  shape  (see  Figure 
12)  it  should  not  be  too  surprising  that  erratic  behavior  was  present. 
Stabilizing  the  plasma  with  a  small  axial  magnetic  field  as  demonstrated  by 
MatherAl  may  be  one  promising  approach  to  improve  Salge's  experiment  and 
will  certainly  be  tested  in  stage  4  of  the  proposed  research  schedule 
(during  the  third  year,  see  Figure  2). 

Another  approach  to  firing  control  might  be  to  allow  the  gas  to 
completely  recombine  between  shots,  h  plasma  focus  has  been  operated  at  a 
rate  of  1  Hz  for  150  shots. ^2  Only  two  shots  in  this  sequence  had  no 
neutron  yield  and  the  rest  were  within  50%  of  5 *  10^  neutrons/shot. 

Since  strong  neutron  yield  is  associated  with  sharp  current  interruption 
this  latter  result  is  especially  encouraging.  The  upper  limit  to  the 
repetition  rate  in  the  gas  recovery  mode  (besides  being  limited  by  the  power 
supply)  is  the  inverse  of  the  deionization  time.  Since  the  operating 
pressures  and  gap  dimensions  are  roughly  those  of  a  hydrogen 
thyratronA3,  the  deionization  time  would  be  about  5  usee.  Hence  the 
fundamental  limit  on  the  rep.  rate  in  this  mode  would  be  ~lo5  Hz. 

Thus  the  high  voltages  occuring  across  the  DPF  during  Salge's 
experiments  and  the  resulting  erratic  behavior  are  neither  surprising  nor 
fatal  to  the  DPF  opening-switch  concept  although  fairly  sophisticated 
circuitry  may  be  required  for  repetitive  operation  if  the  simple  technique 
of  magnetic  stabilization^  is  unsuccessful.  The  gas -recover y-DPF  mode 
circuitry  is  currently  under  study. 

The  last  two  points  indicate  that  possibly  little  of  the  device  current 
was  carried  by  the  pinch  in  Salge's  experiment  and  hence  the  low  percentage 
interruptions.  Much  stronger  current  interruptions  and  neutron  yields  have 
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been  reported^  for  a  device  of  similar  size  and  carrying  similar 

O 

current  which  had  an  I0  of  4* 10*1  A/s  as  opposed  to  2*1010  A/s 
in  the  Salge  experiments^.  High-power  DPF  experiments  driven  by  a 

O 

magnetic  flux  compressor^  (which  has  a  similar  problem  of  low  I0  as 

O 

does  Salge 1 s  inductive  energy  store)  have  reported^  that  a  minimum  I0 

o 

is  required  for  good  focus  performance  and  I0  minimum  is  dependent  on  filling 
pressure.  This  dependence  ranges  from  1.5*10^  A/s  to  10^  a/s 
for  a  filling  pressure  of  10  to  22  torr  deuterium. 

O 

One  might  expect  an  effect  of  the  I0  on  performance  since  the  time 

O 

the  focus  discharge  dwells  in  a  given  region  depends  on  I0.  Long  dwell 
time  could  cause  that  region  to  heat  and  hence  cause  gas  to  evolve  from  that 
region.  If  enough  gas  evolves  behind  the  current  sheath  as  it  moves  away 
(Figure  A3)  a  second  sheath  may  occur  behind  the  first  sheath  and  hence 
limit  lowering  the  current  in  the  subsequent  pinch  caused  by  the  collapse  of 
the  first  sheath.  This  behavior  has  been  found  to  occur  in  plasma  focus 

°  /\  7 

devices  of  relatively  low  I0.  '. 

The  time,  t,  for  the  current  sheath  to  run  to  a  given  distance  Z  along 
the  central  electrode  can  be  shown  to  be 


7 1/2  i/4  r  o  o  o 

_  1  po  24ff2(b  -a2) 

mr-  irinWar 

o  L  u 


where  p^,  is  the  mass  density  of  the  fill  gas,  and  a  and  b  are  the  radii  of 
the  inner  and  outer  electrodes  respectively.  This  formula  yields  the 
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dependences  on  P0  and  I0  observed  in  experimental  and  gives  order  of 
magnitude  results.  It  can  be  seen  that  the  dwell  time  in  the  vicinity  of 

o  o 

the  breach  (Z=0)  is  inversely  proportional  to  I0*/2.  hence  higher  I0 
should  result  in  lower  initial  dwell  times  and  hence  less  gas  evolution  in 
the  breach  of  the  coaxial  gun  (Figure  A2). 

Parasitic  currents  occuring  in  the  breach  of  the  1MJ  plasma  focus  at 
Frascati^  are  believed  to  be  responsible  for  the  lower  than  expected 
neutron  yields  reported  for  that  device  with  respect  to  those  reported^ 
for  Mather's  DPF  of  similar  energy.  The  Frascati  device  operated  with  an 

O  O 

I0  of  1.3*10*2  a/s  whereas  Mather's  device  had  an  I0  of  3 . 3 *  10^2  a/s. 

An  extreme  example  of  the  process  of  parasitic  currents  occuring  in  the 
breach  of  a  plasma  gun  is  the  erosion  rai 1gunA10-A12#  jn  the  erosion 
rail  gun  the  entire  working  gas  evolves  from  the  insulator  and  electrodes  at 
the  breachA*0.  In  these  experiments  I0  and  peak  current  wereA*2 
3*10*0  A/sec  and  120  kA  respectively  which  are  similar  to  those  of 
Salge.  In  the  erosion  railgun,  greater  than  50%  of  the  device  current 
remains  in  the  breach  regional  and  the  percentage  of  total  system 
energy  dissipated  in  the  breach  was  found  to  decrease^  with  bank 
voltage,  V,  and  to  be  independent^  of  the  number  of  capacitors  (bank 

O 

energy).  Since  I0  =  V/Lx  (where  Lx  is  the  total  external  inductance) 
and  since  Lx  was  dominated  by  the  spark  gapA*2  (hence  independent  of 
number  of  capacitors  and  bank  energy)  this  type  of  dependence  would  be 

O 

expected  if  the  increased  dissipation  were  caused  by  low  I0. 


I 


Thus  there  is  considerable  evidence  that  a  minimum  I0  is  required  to 
keep  parasitic  currents  low  in  the  dense  plasma  focus  and  that  in  Salge's 

O 

experiment  I0  may  have  been  too  low.  Also  large  outward  pointing  electric 
fields  have  been  foundA^  crucial  to  good  current  sheath  formation  and 

O 

these  fields  are  directly  proportional  to  I0.  Only  careful  measurements 
of  the  magnetic  field  distribution  in  the  device  (B  probe  measurements)  can 
answer  this  question  and  such  measurement  were  not  reported  by  Salge*^. 

We  have  proposed  these  types  of  measurements  in  the  research  plan  (Figure  2) 
so  the  results  of  this  research  should  be  less  ambiguous. 

Conclusion 

Thus  while  the  DPF  opening  switch  experiments  performed  by  Salge  were 
pioneering  in  nature,  insufficient  data  is  provided  to  determine  whether  the 
performance  was  ever  optimized  (for  capacitive  drive)  or  whether  parasitic 

O 

currents  dominated  due  to  the  low  I0.  To  use  the  results  of  these 
experiments  as  proof  that  a  DPF  opening  switch  concept  is  impossible  is  to 
greatly  overstate  the  case.  However  one  can  use  the  results  of  Salge's 
experiments  to  suggest  better  modes  of  operation  (e.g.,  the  gas  recovery 
mode  and/or  the  magnetic  field  stabilized  mode)  or  to  suggest  more 
comprehensive  experimental  tests  (e.g.,  B  probe  measurements)  and  this  is 
what  we  plan  to  do  (Figure  2). 
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The  answers  to  the  four  questions  listed  in  Section  Cl  provide  the  goals 
required  for  realizing  the  potential  of  the  Plasma  Focus  Opening  Switch 
(PFOS)  concept.  The  ordering  of  the  goals  is  such  as  to  provide  tests  of 
increasing  difficulty  and  yet  at  the  same  time  opening  the  concept  to  a  wider 
range  of  applications  at  each  step  if  the  result  of  the  previous  goal  is 
positive.  Thus  knowledge  of  the  nature^  of  the  current  interruption  in 
the  plasma  focus  (inductive  or  resistive)  leads  to  strategies  to  improve  PFOS 
performance  in  the  subsequent  tests  of  power  transfer,  control,  and  rep. 
rate.  Subsequently  if  little  power  can  be  transferred  out  of  the  circuit 
then  questions  of  control  and  rep.  rate  are  irrelevant.  Since  each  goal 
depends  on  the  results  of  earlier  goals  most  of  the  discussion  will  be 
centered  on  the  first  goal.  The  others  will  be  discussed  in  terms  of  the 
most  likely  results  of  the  earlier  goals  as  perceived  at  present  to  show  how 
this  research  program  could  lead  to  a  fast  reliable  opening  switch. 

The  Nature  of  the  PFOS  Impedance 

It  is  necessary  to  determine  the  nature  of  the  current  interruption  of 
the  plasma  focus  before  potential  loads  can  be  considered  or  ways  of 
improving  performance  can  be  devised.  For  single  shot  experiments  it  is 
desirable  to  match  source  and  load  impedances  to  insure  the  maximum  transfer 
of  energyB2,B3  an(j  hence  switch  impedances  ranging  from  0.5ft  for  exploding 
wire  experiments^  to  1  or  2ft  for  water  transmission  line^5  experiments 
would  be  desired.  For  a  rapidly  pulsed  system  this  may  not  be  the  optimum 
since  time  to  replenish  the  source  energy  may  seriously  degrade  the  rep. 
rate.  However  if  the  energy  density  of  the  inductive  store  can  be  made  a 
hundred  times  greater  than  capacitive  stores,  interesting  pulsed  power 
sources  could  be  envisioned  that  have  rapid  rep.  rates  and  that  are  not 
prohibitively  large. 
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As  an  example  consider  the  circuit  in  Figure  B1  where  the  opening  switch 
Si  is  replaced  by  a  variable  resistor  which  ramps  up  at  a  constant  rate  R 
starting  when  the  power  pulse  to  the  inductive  load  is  desired  (and  S2  is 
closed).  The  power  P  transferred  from  inductor  l_i  carrying  initial  current 
4  to  the  load  inductor  L2  can  be  written82  aS: 


Ll+L2 


2Rlft(e-8t2-e-26t2) 


3  = 


R(L1+L2) 

l1**2 


Bl) 


where  t  is  time  starting  when  R  is  turned  on.  The  total  energy  transferred 
to  the  load  up  to  time  t  can  be  written 


E ( t )  =  -'-i-  \  L1lf[l-(2e 

(4+l2)2  2  1  1 


-8t  -2Bt 


)] 


B2) 


Taking  the  load  to  be  6  nH  which  is  of  interest  in  high  pulsed  power  plasma 

experiments82  and  R  ~  15 On sec  cons^stent  with  parameters  of  the 

Illinois  plasma  focus  the  two  equations  can  be  used  to  predict  the 
performance  of  a  PFOS.  Two  cases  were  considered:  the  present  experiment 
where  4  -  45  nH  and  4  =  560  kA  and  a  high  energy  inductive  energy 
storage  circuit  where  4  ~  145  nH  and  4  =  10  MA.  The  results  are  shown 
in  Table  BI.  The  data  in  Table  I  indicate  that  a  PFOS  operating  in  an 
inductive  energy  storage  system  could  deliver  power  pulses  that  compare  quite 
favorably  with  the  present  SHIVA  II  implosion  experiments  at  Air  Force 
Weapons  Lab  in  terms  of  peak  power82  and  are  only  a  factor  of  three  lower 
in  terms  of  energy  transferred.  Also  Table  bl  illustrates  that  increasing  R 
shortens  the  energy  transfer  time  and  that  measurable  amounts  of  energy  and 
power  could  be  coupled  out  of  the  Illinois  plasma  focus  experiment  assuming 


R  is  not  effected  by  the  transfer  of  energy.  Finally  the  energy  lost  to  the 
focus  circuit  in  the  switch  (i.e.  J  I2Rs  dt)  is  only  1.45  kJ  whereas 
analysis^  of  the  current  waveform  (Figure  3)  indicates  1/2  of  the 
magnetic  energy  in  the  circuit  or  3.5  kJ  is  lost.  This  could  indicate  that 
the  Illinois  Plasma  Focus  has  a  resistance  even  higher  than  that  used  in  this 
study. 

Thus  the  determination  of  the  nature  of  R  is  a  key  to  the  potential  use 

of  a  plasma  focus  as  an  opening  switch.  For  example  if  R  is  due  to  a 

dissipative  effect  then  methods  to  increase  the  dissipation  rate  such  as 

increasing  the  energy  radiated  by  the  use  of  a  high  Z  working  gas  could  be 

•  • 

examined.  If  the  R  (or  L)  is  due  to  sheath  acceleration  then  methods  of 
increasing  the  sheath  velocity  such  as  the  use  of  hydrogen  at  low  pressures 
as  a  working  gas  would  be  the  local  next  step. 

It  should  be  said  at  this  juncture  that  earlier  plasma  focus  (27  kJ) 
circuit  and  current-sheath  velocity  measurements66*87,  indicated  that  the 
maximum  resistive  component  was  150  mft  whereas  the  inductive  (L*)  component 
was  10  to  20  mft  when  the  center  electrode  is  positive.  However  since  the 
total  current  interrupted66  was  only  10  kA  (out  of  530  kA)  as  opposed  to 
280  kA  (out  of  560  kA)  for  the  Illinois  plasma  focus  (IPF),  similar 
measurements  will  have  to  be  performed  on  IPF  to  be  certain  that  the 
resistive  component  continues  to  dominate.  These  measurements  will  be 
discussed  i'<  detail  in  the  next  section,  but  for  the  basis  of  the  discussion 
of  the  rest  of  this  section  it  will  be  assumed  that  the  resistive 
(dissipative)  component  dominates  the  current  interruption  in  IPF  (and  all 
other  PFOS)  as  well. 


the  Energy  into  a  Load 


While  this  is  an  obvious  next  step  after  determining  the  nature  of  the 
IPF  impedance  it  is  interesting  to  speculate  on  how  the  resistive  component 


might  be  effected  by  withdrawing  current  and  power  from  the  plasma  focus  to 
suggest  possible  ways  to  circumvent  such  behavior. 


As  a  possible  scenario,  assume  the  anomalous  resistance  of  a  plasma 

focus  (typically  100X  Spitzer)  is  caused  by  current  driven  instabilities 

which  are  caused  when  the  electron  drift  velocity,  vq,  (carrying  the 

current)  exceeds  some  typical  plasma  velocity  as  the  electron  thermal 

velocity  vT  or  the  ion  sound  speed  Cs.88  if  current  is  transferred 
1  e 

to  the  load  the  resulting  lowering  of  vq  could  cause  these  instabilities  to 
be  stabilized  (and  damped)  and  hence  cause  the  resistance  of  the  plasma  focus 
to  rapidly  decrease. 

Computer  calcul ations8^  based  on  a  one-dimensional  Lax-Wendroff  two- 


fluid  model  including  anomalous  transport  as  outlined  in  reference  B8 
indicate  that  conditions  should  be  favorable  for  current  driven  instabili¬ 


ties.  These  instabilities  could  increase  the  resistivity  of  the  plasma  focus 
to  over  one  thousand  times  Spitzer  resistivity  during  the  current  collapse 
which  is  the  same  order  as  R  ~  Vp/Ip  estimated  above.  If  this  is  truely 
the  dominant  mechanism  for  enhanced  resistivity  as  the  above  analysis 
indicates  then  tapping  the  power  into  the  load  could  considerably  lower  the 
switch  resistance  just  when  high  resistance  is  desired.  Clearly  some  means 
of  circumventing  these  phenomena  would  have  to  be  sought. 


Another  mechanism  that  would  enhance  resistivity  appears  to  be  occurring 
in  the  plasma  focus810.  This  mechanism  is  magnetically  induced  mass 


motion.  Starting  with  the  MHO  equation  of  motion  and  the  generalized  Ohm's 
Law8**  the  following  relation  can  be  obtained:8^ 
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where  the  closed  loop  moves  with  the  fluid.  V  is  the  fluid  velocity;  is 
the  magnetic  flux  enclosed  by  the  loop;  P,  P,  and  a  are  the  mobility,  mass 
density  and  conductivity  of  the  fluid  respectively,  and  the  rest  of  the 
notation  is  standard.  Since  the  integral  around  the  loop  of  a  gradient  is 
zero  the  last  term  in  equation  3)  can  be  neglected.  This  equation 
illustrates  that  decaying  flux  can  be  converted  into  heat  via  Ohmic  losses  or 
into  mass  motion.  Since  the  latter  also  removes  energy  from  the  circuit  it 
also  represents  an  effective  resistance. 

To  induce  mass  flow,  one  of  the  constituants  of  the  fluid  must  be  free 


to  move  across  the  magnetic  field  lines;  a  figure  of  merit  for  this  can  be 
expressed  in  terms  of  the  Hall  parameter  wcjTj,  where  wcj  is  the 
cyclotron  frequency  for  spec'es  j  and  tj  is  the  collision  frequency  for 
that  species.  For  a  typical  plasma  focus  pinch  toc , <1  for  the  ions 
and  uceTe»l  for  the  electrons  so  the  ions  are  free  to  move  across  the 
field  and  be  accelerated  by  a  decreasing  magnetic  flux  into  the  circuit. 
Since  p/cr  =  -1/ (|  e|  ne)  the  ion  motion  is  away  from  the  center  electrode 
when  it  is  positive  and  toward  the  center  electrode  when  it  is  negative.  A 


solid  center  electrode  would  impede  mass  flow  in  the  latter  case  and  this  may 
explain  why  a  strong  shock  wave  moving  away  from  the  center  electrode  is 
observed  only  when  it  is  positive®^  and  why  the  resistance  of  the  focus 
is  ten  times  higher  with  a  positive  center  electrode  than  with  a  negative 
center  electrode.®7  This  effect  might  also  explain  the  generation  of  1 
kJ  of  fast  ions  moving  away  from  the  center  of  electrode  in  the  Illinois 
plasma  focus®14,B15  wi th  a  positive  center  electrode. 


Magnetically  induced  mass  motion  (MIM),  i.e.  the  conversion  of  magnetic 


energy  to  mass  motion  should  only  be  inhibited  if  the  ion  Hall  parameter 
increases  above  one.  Since 

bt3/2 

03  .T.o<- - - - 

ci  i  n  . 

where  nj ,  T^,  B  are  the  ion  number  density,  the  ion  temperature  and  the 
local  magnetic  field,  ion  heating  or  rarefication  could  stop  this  effect. 
However  tapping  the  current  to  a  load  should  also  reduce  B  so  the  switch 
resistance  may  not  be  effected  by  including  a  load.  This  would  be  a  very 
promising  event  in  view  of  the  results  reported  in  Table  BI. 

Thus  the  plasma  focus  load  experiments  might  yield  both  the  practical 
result  of  coupling  energy  to  a  load  using  a  PFOS  and  also  reveal  the 
resistance  mechanism  which  could  have  important  consequences  beyond  the  PFOS 
concept. 

Firing  Control 

If  magnetically  induced  motion  is  the  cause  of  the  observed  resistivity 
in  the  plasma  focus  this  may  be  enhanced  by  counter  streaming  electron  beams 
(traveling  in  the  opposite  direction  to  the  current  carrying  electrons  in  the 
plasma)  to  enhance  the  flux  annihilation  in  the  circuit.  Since  half  the 
magnetic  flux  being  annihilated  is  that  of  the  circuit  this  would  enhance  the 
resistivity  of  the  main  circuit  and  provide  a  means  of  control  of  the 
switching  process.  Alternatively  two  stream  instabilities  produced  by  such  a 
beam  could  also  enhance  the  resistivity.  Strong  interactions  have  been 
observed  between  electron  beams  and  plasma  focus  pinches®^, B17  ancj  a 
factor  of  five  increase  in  neutron  production  was  observed^  when  a  350  kV, 
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32  kA,  30  nsec  electron  beam  interacted  with  the  focus  in  the  counter  current 
direction.  Since  the  enhanced  ion  mass  flow  predicted  by  equation  B3)  would 
produce  more  fusion  neutrons  by  either  a  beam  target  or  a  thermonuclear 
mechanism,  the  observed  enhancement  in  neutron  production  is  consistent  with 
the  magnetically  induced  motion  model. 

For  these  reasons  electron  beam  plasma  focus  interactions  appear  to  be 
promising  means  of  controlling  the  properties  of  the  plasma  focus  opening 
switch  and  thus  electron  beam  focus  experiments  are  an  important  next  step 
after  the  focus-load  experiments.  Permission  to  use  a  250  kV,  80  kA  and  30 
nsec  electron  beam  has  been  obtained  from  the  Gaseous  Electronics  Lab  at 
UIUC.B18 

Repetitive  Operation 

If  the  electron  beam  plasma  focus  experiments  are  successful  and  the 
magnitude  of  the  interruption  can  be  controlled  by  an  electron  beam,  then 
repetitive  PFOS  experiments  would  be  the  final  step.  This  could  be 
accomplished  by  either  shortening  the  center  electrode  so  multiple  pinches 
could  occur  in  the  first  half  cycle  of  the  plasma  focus  current  or  by  driving 
the  plasma  focus  by  an  inductive  energy  store  as  in  the  case  of  Salge.* 

The  lack  of  firing  control  exhibited  in  the  experiments  by  Salge  may  be 
overcome  by  an  electron  beam  controlled  plasma  focus  discharge.  Other 
strategies  were  outlined  in  Appendix  A.  Hence  the  research  program  proposed 
here  could  lead  to  a  design  for  a  high-power  repetitive  controllable  opening 
switch  if  the  results  of  each  step  of  the  program  are  positive.  Since  the 
goals  of  each  step  are  well  defined  the  results  can  be  readily  evaluated  at 
the  end  of  each  step  as  well. 
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Table  BI 

L1 

45  nH 

145  nH 

145  nH 

145  nH 

h 

560  kA 

10  MA 

10  MA 

10  MA 

• 

R 

10^1 /sec 

10^2 /sec 

3.106o/sec 

107O/sec 

PMAX 

0.94  1010W 

12 

3.4*10  W 

5.9*1012W 

10.7* 1012W 

E 

0.65  kj 

248  kj 

248  kJ 

248  kJ 

90t 

t90f 

125  ns 

130  ns 

75  ns 

41  ns 

R(tgo)t 

0. 125  Q 

0.130 

0.225  0 

0.410 

i/2L i1? 

12.5  kJ 

7.25  MJ 

7.25  MJ 

7.25  MJ 

fl2R  dt 

1.45  kj 

425  kj 

424  kJ 

424  kJ 

tThe  subscript  refers  to  that  quantity  when  90%  of  the  energy,  E,  is 
transferred  to  the  load. 
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Inductive  energy  supply  circuit.  The  energy  is  stored  in  the  inductance 
switch  S-,  closed.  Ideally  the  energy  would  be  transferred  to  load  Lg  by 
simultaneously  opening  S,  and  closing  S„. 
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